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A comprehensive one-dimensional model accounting for the effects of heat and mass transfer, chemical 
kinetics, and drying was developed to describe the torrefaction of a single woody biomass particle. The 
thermochemical sub-models depend only on previously determined or measured characteristics, avoid¬ 
ing the use of fitting or tuning parameters and enabling a rigorous energy balance of the process. More¬ 
over, a high temperature drying sub-model is introduced which overcomes the difficulties associated 
with existing approaches to give physically consistent results, smooth implementation, and numerical 
stability. The particle model was validated against experimental data from the literature for intraparticle 
temperature profiles, particle mass and energy yields over a range of particle sizes and reaction temper¬ 
atures. The modeling results describe well the three distinct stages observed during the torrefaction of 
large particles including the heatup, drying, heat release due to exothermic reactions resulting in thermal 
overshoot, followed by thermal equilibrium where conversion is governed by mass loss kinetics. The non¬ 
linear effects of particle size, temperature, moisture content, and residence time on the mass and energy 
yields are quantified and explained. Larger particles exhibit a significant internal temperature gradient 
and strong temperature overshoot especially at the centerline. The magnitude of the overshoot is a func¬ 
tion of the conductivity, particle size, and average heat release rate. Because of the rise in the reaction 
rate, higher temperatures increase the sensitivity of the process to particle size. Due to the dependence 
of drying rate on heat transfer limitations, the sensitivity of torrefaction to initial moisture content 
increases strongly with particle size. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Because of increasing concerns over rising greenhouse gas 
emissions, biomass has received interest as a renewable feedstock 
for a variety of thermochemical processes including combustion, 
gasification, and pyrolysis. However raw biomass displays several 
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unfavorable properties related to its low bulk energy density, high 
moisture content and uptake rate, and exorbitant size reduction 
energy penalty. Torrefaction has been proposed and studied as 
thermal pretreatment approach to improve these characteristics 
[1,2], Torrefaction is a mild pyrolysis process occurring between 
200 and 300 °C over a residence time between several minutes to 
about an hour which results in partial devolatilization (0-60% by 
weight) of the original dry feedstock [3], During torrefaction, the 
release of the volatiles combined with other changes to the physi¬ 
cal and chemical structure result in a solid product with greater 
specific energy density, increased resistance to fungal degradation 
[4], and reduced grinding energy requirement [5], 
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Nomenclature 



A 

pre-exponential factor (s 1 ) 

Greek letters 

B 

permeability (m 2 ) 

A H evap 

enthalpy of evaporation (J kg -1 ) 

Cp 

specific heat (J kg 1 I< ') 

8 

porosity () 

d 

pore diameter (m) 

rienergy 

energy yield () 

D 

Diffusivity (m 2 s ') 

P 

viscosity (Pa s) 

E 

activation energy (J moD 1 ) 

P 

density (kg m 3 ) 

h c 

convective heat transfer coefficient (W m 2 K ’) 

CO 

emissivity coefficient ( ) 

hr 

radiative heat transfer coefficient (W m 2 I< ’) 

a 

Stefan-Boltzmann constant (W m 2 K 4 ) 

H 

total enthalpy (J kg -1 ) 



H° f 

standard heat of formation (J kg -1 ) 

Dimensionless numbers 

HHV 

higher heating value 0 kg -1 ) 

Bi 

Biot number 

k 

thermal conductivity (Wm -1 K ’) 

Nu 

Nusselt number 

k 

reaction rate (s ’) 

Re 

Reynolds number 

L 

particle length scale (m) 

Cr 

Grashof number 

m 

particle mass (kg) 

Pe 

Peclet number 

MW 

atomic weight (kg mol 1 ) 



n 

shape factor () 

Abbreviations 

MC d 

moisture content (% dry basis) 

d.b. 

dry basis 

P 

pressure (Pa) 

w.b. 

wet basis 

Qr 

heat of reaction rate (W m 3 ) 

radial (sphere and cylinder), axial (slab) distance (m) 



R 

R g 

particle size (radius) (m) 
universal gas constant (J mol 1 K ’) 

suDscnpzs 

0 initial condition 

T 

T b 

temperature (K or C) 
vaporization temperature (K) 
time (s) 

j 

i 

hnal condition 

pseudo-component (i.e. A, B, C, 171,172, S) 
chemical specie (i.e. C, H, O, N, S) 

v 

velocity (ms 4 ) 

p 

particle 

Y 

mass fraction (kg/kg) 

s 

gas 


fsp 

fiber saturation point 


Torrefaction is a complex process involving coupled intraparti¬ 
cle physical and chemical phenomena including external heat 
transfer, intraparticle heat transfer, thermochemical decomposi¬ 
tion, convection and diffusion of volatile products. For sufficiently 
small particles (<2 mm), torrefaction is essentially free of heat and 
mass transfer limitations [6,7], The solid mass loss kinetics of tor- 
refaction in this regime has been studied for a variety of woody 
feedstocks including willow [6], beech [8] and softwoods (pine, 
fur, spruce, and pine bark) [9], The torrefaction of the individual 
lignocellulose components (xylan, cellulose, and lignin) has also 
been examined [10], 

Because of the high grinding energy requirement of raw bio¬ 
mass [5], industrial scale torrefaction will preferably be performed 
with large (>2 mm) particles. In this case, significantly more com¬ 
plex interactions occur. Heat transfer limitations cause a tempera¬ 
ture gradient across the particle resulting in chemical reactions 
proceeding non-uniformly. The non-uniform release of heat due 
to exothermic decomposition reactions can further exacerbate 
the heterogeneity of the process and even lead to thermal runaway 
[8,11 ]. In order to examine these phenomena, robust coupled tran¬ 
sient models describing the transport, kinetics, and thermochemis¬ 
try are required. Numerous ID and 2D physico-chemical models 
exist in the literature that describe single particle pyrolysis at high 
temperatures (400-700 °C) [12-15], However, pyrolysis at these 
temperatures results in partial to complete devolatilization (70- 
90 wt%) and different thermal degradation pathways are expected 
compared to torrefaction. 

Fewer models have considered coupled heat/mass transfer and 
kinetics at low temperatures (<300 °C) consistent with torrefac¬ 
tion. Turner and coworkers [11,16] developed a 2-D computational 
heat and mass transfer model for the torrefaction of beech (Fagus 
sylvatica ) wood and compared the predicted temperature profiles 


to experimental results using large (5 x 15 x 25 cm) beech boards. 
Their assumed kinetic mechanism described wood thermal 
decomposition as a linear superposition of the decomposition 
hemicelluloses, cellulose, and lignin components. Turner noted 
that the reaction enthalpy term has a very important effect on 
the temperature profiles [16], however due to the lack of available 
parameters in the literature, the enthalpy of reaction of the hemi- 
cellulose decomposition was used as a tuning parameter. Kadem 
and coworkers [17] developed a transient 3-D computational 
model of the heat treatment (150-250 °C) of wood, focused on heat 
and moisture transfer (drying), and neglected thermal decomposi¬ 
tion kinetics. Patuzzi and coworkers developed and validated a 
model describing the heat, mass, and momentum transfer during 
the torrefaction of common reed (Phragmites australis ) in a custom 
bench-scale reactor [18], Parameters for a two stage solid mass- 
loss kinetics mechanism were separately fitted, and the global heat 
of reaction (including both drying and torrefaction) was treated as 
a calibration parameter. Ratte et al. [19,18] developed a particle- 
scale model coupled with a reactor model to describe their labora¬ 
tory pilot (3-8 kg/h) moving bed reactor system [19,20], They 
adopted a global kinetics decomposition model while the thermo¬ 
chemical model postulated a fixed composition of solid and vola¬ 
tile products enabling them to explicitly model their enthalpies. 
Although the particle sub-model [19] was not compared explicitly 
with experimental particle temperature profiles, the reactor model 
was validated against measured gas temperature profiles [20], 
Recently, Basu and coworkers [21] compared a ID particle scale 
model with experimental measurements of solid mass loss curves 
and centerline temperature profiles from cylindrical poplar 
dowels. While the solid mass loss kinetics of the first stage were 
separately fitted under TGA conditions, multiple kinetic and 
thermochemical parameters were tuned during the development 
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of the particle scale model including the heat of reaction. In 
reviewing the previous modeling efforts, no models exist which 
simultaneously (i) show detailed validation with intraparticle tem¬ 
peratures, (ii) are free of tunable parameters, and (iii) give informa¬ 
tion on the energy yield of the solid product. 

In this paper, we formulate a particle scale model that couples 
drying, torrefaction kinetics, and relevant transport processes 
which overcomes these aforementioned deficiencies. Ultimately, 
the model is then used to examine the effects of particle size, tem¬ 
perature, and moisture content on the overall rate and uniformity 
of the process. We use previously determined thermochemical and 
kinetic data and validate the model against published experimen¬ 
tal data. In such a way, the model is formulated without the use of 
adjustable or fitting parameters. In the first part of this paper, a 
characteristic time scale analysis is applied to identify the domi¬ 
nant physical and chemical phenomena occurring during torrefac¬ 
tion of large particles. Next, the mathematical formulation of the 
coupled transport and kinetics model is presented along with the 
selected kinetics and thermochemistry model. In the following sec¬ 
tion, the results are validated against the experimental measure¬ 
ments of Van der Stelt [8] and Basu et al. [22], The discussion 
highlights the models ability to capture important behaviors 
occurring during single particle torrefaction including the effects 
of particle size, moisture content, temperature, and exothermic 
overshoot. 


2. Physical and chemical processes during torrefaction 

Biomass particles undergoing torrefaction experience a com¬ 
plex set of interacting phenomena including mass and heat trans¬ 
fer as well as drying and chemical devolatilization. It is important 
to judiciously identify the most relevant processes in order to for¬ 
mulate representative models. 

2.1. Drying 

Drying is one of the first processes to occur during the torrefac¬ 
tion of a biomass particle and constitutes a significant thermal load 
during the process. After an initial heating period, drying com¬ 
mences once the local moisture boiling temperature is reached, 
which is a function of moisture content and ranges between 100 
and 135 °C [23,24], At a particular location in a particle, drying will 
be completed before torrefaction initiates because temperatures 
greater than 200 °C are required for the latter. However, if heat 
transfer limitations and intraparticle temperature gradients are 
present, these processes may overlap whereby torrefaction occurs 
on the hotter particle exterior while the cooler interior is still heat¬ 
ing up or drying. 

Moisture present in fresh wood varies between 15 and 60% by 
mass on a wet basis (w.b.) depending on the type, location and 
method of storage [25,26], Even after the application of a thermal 
drying process, the equilibrium moisture content will still remain 
above 6-10% w.b. depending on the relative humidity and temper¬ 
ature [27-29], Water exists in three forms within biomass includ¬ 
ing (i) vapor, (ii) free (or capillary) water, and (iii) bound (or 
hygroscopic) water [28], Bound water is adsorbed through hydro¬ 
gen bonding [30] to the cell wall constituents (hemicelluloses, cel¬ 
lulose). Once the adsorption sites are saturated with bound 
moisture, the fiber saturation point is reached (typically 23% w.b. 
or 30% d.b.) [25], Above this moisture content, water fills the voids 
or pores of the biomass and is referred to as free water. 

During torrefaction, free and bound moisture can be completely 
evaporated and resulting in a product with as low as 0-2% mois¬ 
ture [31], Additionally, the decomposition of hydroxyl groups in 
the cell wall components resulting in a more hydrophobic product 


[32,33], While the effect of initial moisture content on the torrefac¬ 
tion process has yet to be studied experimentally, drying has pre¬ 
viously been shown to affect pyrolysis times [28], Because of the 
large latent heat of vaporization of liquid moisture, the thermal 
load associated with drying - even for a relatively dry particle of 
6% w.b. - will be significant (-30%) compared to the sensible heat 
required to raise the temperature of the particle from ambient to 
reaction temperatures (200-300 °C). 

The dominant drying mechanisms depend on temperature and 
can include bulk convective and diffusive flow of free and bound 
moisture. Under high temperature drying conditions (>100 °C) 
the diffusive flow of bound moisture and vapors can be considered 
negligible compared to the convective flow of vapor [34,23], 

2.2. Characteristic time scale analysis 

A characteristic time scale analysis similar to that performed by 
Chan et al. [12] aids in identifying the dominant processes occur¬ 
ring during the drying and torrefaction of a single particle and is 
summarized in Table 1. Although, the time scales do not represent 
the magnitude of the driving forces behind such processes (i.e. 
temperature gradient) they allow the identification of the rate con¬ 
trolling mechanisms. The analysis provides justification for the 
approximations and assumptions made in the mathematical model 
presented next. Values for wood and gas physical properties - 
including conductivity (fc), specific heat (c p ), and density (p) - were 
taken from [35], The value for the diffusivity (D), was estimated by 
the correlation in [30] while permeability (B 0 ). and gas viscosity (p) 
were taken from [14], The solid mass loss reaction rate was evalu¬ 
ated with the willow mass loss kinetics at 250 and 300 °C [3] and 
the drying rate was evaluated at 100°C [36], The length scales 
described in Table 1 include the micro scale, which corresponds 
to the pore diameters, and the macro scale - reflecting the maxi¬ 
mum size of feed particles. 

Several important observations can be drawn from Table 1. 
Mass transfer, especially pressure-driven transfer, occurs at time 
scales much faster than either heat transfer or reactions. At the 
micro scale of the wood pores (tracheids), convective heat transfer 
is rapid compared to kinetics, and the heat capacity of gas is low 
compared to that of the solid. As a result, gas leaving the solid rap¬ 
idly heats up to the temperature of the wood/char it flows through. 
If a positive temperature gradient exists in the particle (i.e. the 
exterior of the particle is hotter than the interior of the particle), 
these outward flowing volatiles can cool the outer hotter layers 
as they leave. This effect is often neglected for simplicity 
[21,35,37] but has been shown to be significant, especially during 
drying [34,38], The relative importance of the convective term to 
the internal conduction term can be quantified by the relevant Pec- 
let number. Under conditions of active drying or devolatilization, 
the Peclet number has an order of magnitude of 1, confirming that 
the heat transfer associated with the flowing volatiles should be 
included. 

Pe = Lp g cp g u/k som . (1) 

As expected, torrefaction kinetics have a strong temperature 
dependence occurring almost an order of magnitude faster at 
300 °C than 250 °C. Drying occurs relatively fast compared to heat 
transfer at the particle scale. Considering the characteristic length 
scale of relatively large particles (L = 0.01 m), mass loss kinetics, 
external heat transfer (convection), and conduction have time 
scales of similar order of magnitude. Equivalently, the Biot number, 
which is the ratio of the characteristic time for internal heat trans¬ 
fer divided by the time for external heat transfer is approximately 
1. Therefore no rate-limiting simplifications can be made and the 
model must account for these transport mechanisms. 
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Characteristic times of physical and chemical processes occurring during torrefaction. 


Transport process or reaction 

Characteristic time (s) 

T= 250-300 °C 


Microscale I = 10 pm 

Macroscale I = 1 cm 

Diffusion, i 2 /D eJf 

10 4 

10 2 

Intraparticle fluid flow, p! 2 /ApB 0 

10 5 

10’ 

Convective/radiative heat transfer, (pc p ) s I/h 

10° 

10 3 

Conduction heat transfer, (pc p )J?lk 

10 3 

10 3 

Torrefaction mass loss kinetics, 1 /k, 

10 2 -10 3 

10 2 -10 3 

Drying kinetics 1 /k d 

10 2 

10 2 

Deff*= 10- 6 m 2 s- 1 

Bo ss 10~ 14 m 2 


(pCp) s m 10 6 J m 3 K 1 

hK 10Wm 2 r' 


(pc p ) g s3l0 4 Jm- 3 K-' 

kralO 1 Wm ’ K 1 


p « 11T 5 Pa s 

k, ss 239, 1300 s 1 (at 300, 250 °C) 


Ap S310 4 Pa 

k d ss 440 s 1 



A summary of the conclusions derived from Table 1 : 

(1) Mass transfer via diffusion is much slower than that by 
hydrodynamic (pressure driven) flow from drying or 
devolatilization. 

(2) Mass transport by convection is rapid and occurs much fas¬ 
ter than heat transport, drying, and chemical devolatiliza¬ 
tion processes. 

(3) The thermal heat capacity of the gas ( pc p ) g is much less than 
that of the solid {pc p ) soUd . As a result local thermal equilib¬ 
rium exists between the volatiles and the solid. Internal con¬ 
vective heat transfer is significant under conditions of active 
devolatilization. 

(4) Significant heat transfer limitations exist at the particle 
scale, which are neither controlled entirely by external or 
internal heat transfer. 

3. Mathematical model 

The formulation of the numerical model is based on the follow¬ 
ing assumptions: 

(1) Particles are represented in a one-dimensional, time-depen¬ 
dent domain. A shape factor enables the description of slab, 
cylinder, or spherical geometries. 

(2) Total volume of the particle remains unchanged during tor- 
refaction. In other words, no overall structural changes (i.e. 
shrinkage, breakage) occurs. 

(3) Outflow of volatiles is instantaneous such that the pressure 
of the particle remains unchanged during drying and 
torrefaction. 

The physical processes described by the model include: 

(1) Heat transfer by radiation, convection, and conduction to the 
exposed surfaces of the particle and within the solid. 

(2) Convective mass transport of the volatiles through the 
particle. 

(3) Variable thermal and physical properties (moisture content, 
conductivity, density, heat capacity) with position and time. 

Biomass is inherently anisotropic with properties like thermal 
conductivity and permeability being especially dependent on the 
grain orientation [30], However, the assumption of one-dimen¬ 
sionality along with an effective set of physical properties is valid 
and commonly made in the literature [14,23,37,39-41] during 
the uniform heating of cylindrical biomass particles of sufficiently 
long aspect ratio. Acceptable agreement between ID models and 


experimental data with cylindrical particle aspect ratios as low 
as 2 for Galgano and Di Blasi [40], 5-8 for Lu et al. [14] and 2.7 
for Pyle and Zaror [37], Alves and Figueiredo [23] found that 
beyond an aspect ratio of 3, the pyrolysis time was hardly affected. 

3.1. Drying 

Three existing modeling approaches have been applied for esti¬ 
mating the evaporation rate of large biomass particles undergoing 
high temperature drying, pyrolysis, or combustion. These include 
the (i) first order kinetic evaporation model [12,34,36,42-44], (ii) 
equilibrium model [30,45-47], and (iii) a constant temperature 
model (also known as the heat sink model) [23,40,41,48-50], The 
first order kinetic approach is a simplified approach which treats 
the evaporation rate as a thermally activated process described 
with a first order Arrhenius rate equation. The activation energy 
and pre-exponential factor are chosen to produce an evaporation 
rate which increases rapidly at the boiling temperature [34,51], 
The advantages of this method are that it is implemented easily 
as an additional chemical reaction in particle models and results 
in smooth intraparticle gradients [52], The disadvantages are that 
some drying occurs before the boiling temperature is actually 
reached, and the kinetic parameters in the literature show a wide 
scatter [53], 

The equilibrium approach accounts for heat and mass transfer 
limitations and assumes that the bound and free liquid moisture 
exist in equilibrium with the local gas phase. The partial pressure 
of water vapor is fixed by the local saturation pressure. The disad¬ 
vantage of the equilibrium approach is that it is known to result in 
depressions in pressure in locations ahead of the drying front and 
requires simultaneous solution of algebraic and differential equa¬ 
tions [52], 

The heat sink approach assumes that drying is controlled 
entirely by heat transfer: when a moist particle layer has reached 
the boiling temperature (e.g. 100 °C), all energy flowing into that 
layer is consumed entirely by evaporation of water until drying 
ceases. It is implemented in models in one of two ways. Either 
the drying front is assumed to be infinitely thin and the particle 
is divided into moist and dry sections or the governing equations 
must include a conditional statement. The former method, though 
easily implemented, is not necessarily valid in cases where the dry¬ 
ing front has a thickness comparable to the particle [44], In the 
case of torrefaction where the Biot number is of order 1, the drying 
front will not be thin. The latter method requires special numerical 
treatment and is not computationally efficient due to the sharp dis¬ 
continuities it introduces into the governing equation [52,53], 

In the current model, an approach which combines aspects of 
the existing methods is adopted where the drying rate is evaluated 
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by a first order Arrhenius equation but has a conditional depen¬ 
dence on temperature (see Eq. (6)). The drying rate is set to 0 for 
temperatures below the boiling temperature. As done by Alves 
and Figueiredo [23] and Bilbao et al. [49] who adopted the heat 
sink approach, the boiling temperature is set by the local moisture 
content and is fitted to data from [24], This approach solves two 
major drawbacks of the first order Arrhenius approach. The drying 
profile is more physical because drying does not occur until the 
boiling temperature is actually reached. Also, it is less sensitive 
to the assumed kinetic parameters because the boiling tempera¬ 
ture is specified. At the same time, it is extremely easy to imple¬ 
ment as an additional chemical reaction without modifying 
governing equations and does not entail the solution of additional 
algebraic equations. 

3.2. Solid mass loss kinetics 

A variety of models have been proposed and verified to describe 
the mass loss kinetics of during torrefaction. Approaches to model¬ 
ing the mass loss kinetics can be divided into two main groups. 
One approach is to describe the wood/biomass as three indepen¬ 
dently reacting constituents (hemicelluloses, cellulose, and lignin) 
[16,20]. This approach is adopted when ease of application to dif¬ 
ferent feedstocks is the goal. Nocquet et al. recently assessed the 
additive approach and found that it reproduces solid mass loss 
curves up to temperatures of 250 °C; however, above this temper¬ 
ature, significant interactions between the constituents were 
apparent - resulting in a large discrepancy [54], The second more 
common approach is a lumped approach where the wood is mod¬ 
eled as one or multiple solid pseudophases reacting to volatiles and 
char. Several authors [8,9,55] have found that wood torrefaction 
could be modeled with two independently reacting components 
forming two volatile phases and two char phases. The mechanism 
adopted in the present study was originally proposed by Di Blasi 
and Lanzetta [56] to describe hemicellulose (xylan) mass loss 
kinetics but was later fitted to satisfactorily describe willow [3], 
common reed [18], eucalyptus, and spruce [55] for temperatures 
up to 300 °C. In this mechanism, the torrefaction products are 
lumped into five pseudo-components and decompose according 
to first order kinetics in two steps with the first step being much 
faster than the second as follows: 



where the solid phase pseudo components (A, B, C ) represent raw 
biomass (A), an intermediate solid (B), and char (C), volatile prod¬ 
ucts are represented by VI and V2, and k u k 2 , kvi. k v 2 represent 
Arrhenius kinetic parameters. 

k i =A i exp(-E i /R g T). (2) 

The pre-exponential factor has units of s \ the activation 
energy has units of J mol \ R g is the universal gas constant in 
J mol 1 1< '■ and T is temperature in K. See Table 2 for the parame¬ 
ters. It was shown that the composition of the lumped volatile 
pseudo-phases (1/1 and V2) were individually well described as 
mixture of nine identifiable chemicals: acetic acid, water, formic 
acid, methanol, lactic acid, furfural, hydroxyacetone, carbon diox¬ 
ide and carbon monoxide [57], See Table 2 in [57] for the details 
of these compositions. The advantage of this kinetic mechanism 
and set of parameters is that all pseudo-components (A, B, C, 1/1, 
1/2) have well-defined chemical compositions thus enabling esti¬ 
mation of their thermodynamic properties, namely their heat of 
formation and specific heat capacity [58], Thus, the elemental 


balance enables the estimation of the solid product composition 
and energy yield [59] not considered in previous single particle 
models. 


3.2.1. Conservation of species/mass 

The equations describing the evolution of the solid pseudo com¬ 
ponents (A, B, C) and liquid moisture ( MC ) are: 


d ^ = -(ki+kvi)p A , 

(3) 

d ^ = kip A -(k 2 + k V2 )p B , 

(4) 

d 3F = k2PB ' 

(5) 

d -^ = -k iPMC , if T > T b , 

(6) 


where p; is the apparent density of pseudo component (i = A, B, C, 
IVfC) in (kg m 3 ) and T b is the vaporization temperature. The condi¬ 
tionality of the evaporation rate on temperature prevents drying 
from occurring before the boiling temperature has been reached. 
The dependence of vaporization temperature on the dry basis mois¬ 
ture content is fitted to the data from Alverez Noves and Fernandez- 
Golfin Seco [24], The local dry basis moisture fraction is given by, 


MCd = Pmc 
P s 


Am 

Pa + Pb + Pc ' 


(7) 


The total solid phase (dry) is modeled as a mixture of the three 
solid pseudo components, 


dp s dp A dp B dp c 
dt dt + dt + dt 


-(fcvi p A + k V2 p B ). 


(8) 


The elemental composition (i.e. the mass fraction of carbon, 
hydrogen, oxygen, nitrogen and ash) of the solid phase on a dry 
basis is given by, 

d ^M« -kv, p A Yj,v 1 - k V2 p B Y m , (9) 


where Y, ,- is the mass fraction of elemental species “j” (j = C, H, O, N, 
Ash) contained in component “i” (e.g. S, 1/1, V2). See Table 2 of [57] 
for these mass fractions. The conservation of mass on the total gas 
phase is given by, 


d(spg) 

§1 


r" dr 


{f'PgV) m 


dPs 

dt 


9p M c 

dt 


= kvi p A + k V2 p B + k d p MC , 


(10) 


where n is a shape factor describing the particle n = 1, 2,3 for a slab, 
cylindrical, and spherical geometries. As was previously done in 
several other single particle pyrolysis models [12,23,60], the mass 
flux of volatiles (kg/m 2 s) can be evaluated from this equation 
because of the assumption of fast volatile release which results in 
a negligible accumulation term: 

P g v = ^ r n (kvtP A + k V2 p B + k d p MC )dr. (11) 


3.3. Conservation of energy 


The conservation of energy for the solid phase accounts for 
accumulation, conduction convection, and reaction, 


gj (Pa c pA + PsCp,B + p c c PtC + sp g cp g ) 


dT (dk kn 
dr\dr + T 


(12) 


where the heat of reaction term, 
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Assumed thermodynamic and physical properties for single particle torrefaction simulation. 


Property 

Correlation/value 

Source 

Dry wood density 

p A o = ps 0 = 700, 500 kg nr 3 (Beech, Poplar) 

[30,63] 

Initial wood composition (dry basis) 

Y JS0 = 0.4692, 0.0585, 0.4662, 0.0001, 0.006 (Beech) for j = C, H, O, N, Ash 

[8] 


Y js o = 0.4845, 0.0585, 0.4369, 0.0047, 0.0143 (Poplar) 

[22] 

Solid specific heat capacity 

cp, =* [e”(^)‘ 2 + 2 ,jkg-’ K’ 

[65] 

Solid average molecular weight 

MW, = |jjg j^lJ i = A, B, Cj = C, H, 0, N kg moV’ 

HHVt = 1000[351.69(Y C ) + 1162.46(Y„) - 110.95(Y o ) + 104.67(Y S ) + 62.8(Y„)] J kg-' 

[65] 

Solid higher heating value 

[66] 

Solid standard heat of formation 

% = HHV t + + Jfc H;, COj t = A, B, C. S j kg-’ 

[58] 

VI specific heat capacity 

cp^ = 1162 + 0.961 TJ kg-’ K ’ 

[58] 

V2 specific heat capacity 

cpv2 = 601 + 2.262TJ kg-' K ’ 

[58] 

VI standard heat of formation 

H° fv , = -10,345,151 Jkg-’ 

[58] 

V2 standard heat of formation 

H) yl = -7,078,295 Jkg-' 

[58] 

Heat of vaporization 

A H° tmp = 2,260,000 J kg-’ for MC d > M fsp 

[40] 


A H° e va p = 10 6 [3.348 - 13.085(MC d ) + 60.262(MC d ) 2 - 95.778(MC d ) 3 ] J kg ' for MC d < M fip 

[49] 

Wood thermal conductivity 

k A = 0.209, 0.11 W m ’ K 1 (Beech, Poplar) 

[14,30] 

Torrefied phase (B) thermal conductivity 

k B = (ka + kc)/2Wm-’K-’ 

(Estimated) 

Char (C) thermal conductivity 

k c = 0.071 W m-’ K ’ (Beech, Poplar) 

[14,30] 

Moisture thermal conductivity 

k MC = (0.406MC d )pso/1000 W m“’ K ' 

[44] 

Effective thermal conductivity 

ke ff = ( p A k A + p B k B + pckc)lpAo + kMc + sk^ + 13.5<rdr7co W m-’ K 1 

[12] 

Gas thermal conductivity 

k gos = 25.77*10- 3 Wm’K-’ 

[38] 

Emissivity coefficient 

co = 0.95 

[35] 

Porosity 

e = 1 - (1 - boXPa + p B + Pc)/Pao 

[15,35] 

Pore diameter 

d = 2* 10 5 m 

[35] 

Stefan-Boltzmann constant 

a = 5.67 * 10 8 W m 2 K 4 


Kinetic rate parameters 

k\ (A -> B) 

2.48 * 10 4 exp(-75,976/RT) s ' 

[36] 

k„ (A -VI) 

3.23 * 10 7 exp(—114,214/RT) s -1 


kz(B^C) 

1.1 * 10’° exp(—151,711/RT) s“’ 


kv2 (B - V2) 

1.59 * 10 lo exp(—151,711/RT) s ' 


Drying rate parameters 

k d (MQ lq -» MC vap ) 

4.5 * 10 3 exp(—45,000/RT) s ' 

[6] 

Vaporization temperature 

T„ = 73.204 exp (26.451 MC d ) + 373 K 

Fitted from [24] 


a-i E If < 13 > 

includes the heat of vaporization AH evap of liquid moisture inj kg -1 . 
Due to the heat of sorption of bound moisture, the heat of vaporiza¬ 
tion term shows a dependence on dry basis moisture content (see 
Table 2) [30,49], The total enthalpy at temperature T is given by 
the summation of the formation enthalpy and sensible enthalpy, 

Hi(T) = HJf + j\ P) i(T)dT, (14) 

where ; is the standard heat of formation of species i inj/kg and c Pii 
is the specific heat capacity of component i inj kg -1 1< \ T is the reac¬ 
tion temperature in Kelvin, and T 0 is the standard temperature 
(298.151<). The total gas phase and VI and V2 are assumed to behave 
as an ideal gas mixture of their nine constituent species. Thermody¬ 
namic properties of the solid phases (A, B, C, S) and VI and V2 are sum¬ 
marized in Table 2. For further information on the thermodynamic 
property estimation of the solid and volatile products see [58], The 
effective conductivity of the solid matrix described in Table 2 
includes a contribution from radiation heat transfer through the pore 
spaces [12] as well as a contribution from the moisture [44,49], 

3.4. Initial and boundary conditions 

The initial conditions for t = 0: 

T = T 0 , p A = p A o, p B = Pc = 0, PJts = PsqYjso- (15) 

The initial chemical composition of the particle (Y c $ o. Yfcso 
• • ■ Ytehs o) depends on the feedstock and is summarized in Table 2. 


The thermal boundary conditions for the particle are shown sche¬ 
matically in Fig. 1. The Neumann boundary condition for t> 0 at 
the centerline r = 0 reflects symmetry, 


The Neumann boundary conditions at the surface of the particle 
r = R includes both convective and radiative heat transfer, 

kg=h c (T-T r ) + OK7(T 4 -7t), (17) 

where h c is the convective heat transfer coefficient in W m 2 K 1 
and T r is the reactor temperature in K. The temperature profile of 
the reactor is not coupled with the particle and is prescribed 
according the conditions of the particular case of interest. In the 
case of Van der Stelt’s experiments, the reactor temperature 
increased at a heating rate between 8 and 16 K/min until reaching 

Particle surface boundary conditions 


= h c (T - T r ) + a>o(T* - T r 4 ) 



T(r, t) = Particle temperature (K) 
T r (t) = Reactor temperature (K) 


Fig. 1 . Particle thermal boundary conditions. 
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the final temperature after which it remained constant [8], In Basu’s 
experiments [22] the particles were dropped into an already heated 
reactor. For the axial flow of air over stationary cylinders the 
convective heat transfer coefficient can be estimated using the 
correlation by Nowak and Stachel [61], 

Nu=^= 0.87 Re p 052 + 1.75Re° 21 + 0.5(0.5Re^ 75 + 8*10 9 Gr* 24 )°° 8 , 

(18) 

where Re p and Gr p represents the Reynolds and Grashof num¬ 
bers, respectively, for the particle. This correlation is valid for 
2 < Re p < 1000 and 0<Gr p <10 7 . Other correlations exist in the 
literature for other ID geometries (i.e. sphere, slab) [14], Substi¬ 
tuting in quantities for the flow and temperature characteristics 
of the experimental reactors under consideration, the convective 
external heat transfer coefficient for cylindrical particles 
between 4.76 mm and 25.4 mm diameter ranges between 20 
and 9 (W m 2 K ’), respectively, decreasing with increasing par¬ 
ticle size. The effective radiative heat transfer coefficient is of 
similar order of magnitude ranging between 20 and 40 
(W m 2 K 1 ), depending strongly on reactor and particle surface 
temperature, thus justifying its inclusion in the boundary 
conditions. 


3.5. Numerical method 

The system of partial differential equations is solved using the 
method of lines approach where finite differences are used to 
approximate the spatial derivatives. The resulting system of ordin¬ 
ary differential equations is integrated using odel5s in MATLAB. 
Treatment is required for the Heaviside step function introduced 
by the conditionality of evaporation rate on temperature. An error 
function was chosen to approximate the step function which was 
found to provide good accuracy while having a minimal impact 
on the integration time. 


4. Results and discussion 

4.1. Comparison with the experiments of Van der Stelt 

4.1.1. Description of experiments 

Van der Stelt performed torrefaction experiments using large 
cylindrical beech and willow particles with diameters between 1 
and 2.8 cm and lengths of 10 cm. The particles were placed in a 
vertical split tube oven which was electrically heated. Argon was 
used as an inert gas to remove the produced volatiles from the 
reactor. Five thermocouples embedded inside the particle continu¬ 
ously monitored the temperature at radial positions (r = 0, 0.3, 0.7, 
1, and 1.3 cm). The particle and reactor - initially at ambient tem¬ 
perature - were heated at a rate of approximately 10 °C/min to the 
final reactor temperature, which was varied between 230 and 
300 °C. The mass and/or energy yield of the particle was not 
reported. The initial moisture content of the particles was reported 
as approximately 6%. 

4.1.2. Intraparticle temperature profiles 

In Fig. 2a and c the experimentally measured temperature pro¬ 
files of the beech particle at two radial positions (r = 0 and 
r= 0.013 cm) are compared with the model predictions for two 
reactor temperatures, 235 and 289 °C, respectively. Hereafter these 
positions are referred to as the centerline (r = 0) and the near sur¬ 
face (r= 0.013 cm) positions. In both cases, the temperature of the 
reactor does not heat up to the final temperature in a perfectly lin¬ 
ear fashion and for modeling the boundary conditions, the reactor 
profile is approximated by two separate linear heating periods. In 
Fig. 2a where Tfi na i = 235 °C, the time for the center of the particle 
to reach 100 °C is 16.5 min, which is the same as the value reported 
in the experiment demonstrating good agreement on the initial 
heating period and drying dynamics. During the drying period, 
the temperature of the centerline of the particle is significantly 
higher than 100 °C (the boiling temperature for liquid water at 
atmospheric pressure). In fact, the moisture at the centerline pla¬ 
teaus at ~114 °C. This demonstrates that the dependence of boiling 
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temperature on moisture content is crucial to representing the 
location of the drying plateau. Drying models which assume a sin¬ 
gle boiling temperature, commonly done when applying the heat¬ 
sink approach, would fail to accurately represent this characteris¬ 
tic. Up until 29 min, the modeled centerline temperature is within 
3 °C of the experimental data. In addition to showing quantitative 
agreement with the Van der Stelt’s results, the drying dynamics 
predicted by model share qualitative similarities with previous 
studies focused on drying of wood cylinders. In Di Blasi’s 
experimental study of moist (50% d.b.) pine cylinder drying at high 
temperatures 200-327 °C [34], the evolution of the intraparticle 
temperature showed a period of initial heatup with low evapora¬ 
tion rates and relatively small particle gradients followed by a dry¬ 
ing period where thermal gradients reached a maximum. After the 
centerline has finished drying, the particle then rapidly approaches 
the final reactor temperature. The same stages are apparent in the 
present results except that after drying has completed, exothermic 
reactions raise the centerline temperature in the model until they 
actually exceed the surface temperature of the particle. 

The temperature profile of the near surface position (r= 0.013) 
predicted by the model is somewhat higher than that experimen¬ 
tally measured profile throughout most of the period, but is within 
the repeatability and uncertainty range for particle temperature 
measurements of this type. While not reported explicitly by van 
der Stelt, in other similar single particle experiments the uncer¬ 
tainty is reported to be at least ±15 °C [13], Particularly, there is 
greater uncertainty near the surface of the particle where large 
temperature gradients exist. The thermocouple - although thin - 
has finite thickness (1 mm) and - besides the possibility that it 
may affect the local effective conductivity - its exact position 
may not be known precisely. 

At the final temperature (at t = 60 min), the centerline is pre¬ 
dicted to be at 241 °C which compares well with the experimen¬ 
tally measured centerline temperature of 240 °C. At 60 min, the 
rate of change of temperature at the centerline and surface of the 
particle is small in both the model and the experiment, demon¬ 
strating that the system has reached a quasi-steady state. At this 
moment, the heat loss from the particle balances any heat released 
from reactions. 

In Fig. 2b for the case of T fina i = 289 °C, the time for the center of 
the particle to reach 100 °C is 13 min for both model and experi¬ 
ment, showing good agreement on the initial heating and drying 
period. A drying temperature plateau at the centerline is well rep¬ 
resented by the model and occurs between 112 and 116 °C. 

The model predicts that the centerline reaches a maximum 
temperature of 329 °C at 38 min which is also in good agreement 
with the experimentally measured maximum temperature of 
335 °C at 39 min. The model predicts that after reaching the peak 
temperature, the centerline temperature decreases at a rate faster 
than the experiment. This is attributable to the endothermic reac¬ 
tions of the second stage which act to cool the particle down. As a 
result at 60 minutes the predicted centerline temperature is 285 °C 
which is 20 degrees lower than the experimentally measured cen¬ 
terline temperature of 306 °C. There are several possible reasons 
for this discrepancy. First, both the solid mass loss kinetics mech¬ 
anism and thermochemical submodel were developed using tem¬ 
perature data below 300 °C, specifically for willow feedstock 
[6,58], This commonly used temperature limit is chosen to avoid 
excessive decomposition of the relatively carbon-rich lignocellu- 
lose components-namely cellulose and lignin [6,57], At tempera¬ 
tures exceeding 300 °C, more complex degradation pathways 
may become significant such as exothermic char-forming reactions 
[62], secondary reactions between the volatiles and the char prod¬ 
ucts [8] and thermal tar cracking. Given the high sensitivity of the 
rates and degradation pathways to temperature, the application of 
these mechanisms 30-40 °C above the original range of fitting 


would strictly be considered extrapolation. Unfortunately, no tor- 
refaction mass loss kinetics mechanisms or thermochemical sub¬ 
models currently have been fitted at higher temperatures (300- 
350 °C), so it may be suggested that future studies may consider 
data in this temperatures range. 

Secondly, the uncertainties of specific modeling inputs may also 
explain some discrepancy during this stage. Namely, in the ther¬ 
mochemistry submodel adopted in this study, the heat of reaction 
of the slower second stage is endothermic but displays more 
uncertainty compared to the first stage which has been positively 
identified as exothermic [58], The uncertainty arises from the var¬ 
iability in the heating values of the solid products as predicted by 
empirical correlations. It was previously suggested that better 
characterization of these products would help to identify whether 
the second stage is either endothermic or exothermic; however, 
the under-prediction of particle temperature by the current model 
suggests the heat of reaction of the second stage is thermally neu¬ 
tral or even slightly exothermic. The sensitivity of the predicted 
temperature profile during this period to the mass loss ldnetics 
and thermochemistry are in agreement with conclusions made 
from previous single-particle models [16], Other physical proper¬ 
ties of the torrefied wood, particularly those which directly affect 
the thermal diffusivity (i.e. density, conductivities, and heat capac¬ 
ities) also display uncertainty but to a lesser degree than the ther¬ 
mochemistry. Whenever the physical property was directly 
measured by the experimentalist, such as wood density, ultimate 
analysis, and initial moisture content, that value was used for the 
respective model simulation (see Table 2). Certain properties, such 
as thermal conductivities and char densities, were not directly 
measured by the experimentalist so in these cases highly cited 
sources from the literature were used to find the appropriate val¬ 
ues and correlations. 

Lastly, another phenomena which could affect the results dur¬ 
ing this stage may be shrinkage effects. While volumetric shrink¬ 
age is negligible during practically relevant torrefaction 
conditions of mass loss less than 20%, at extremely high mass loss 
levels exceeding 50%, the reduction in apparent volume is greater 
than 30% [63], A non-uniform reduction in particle volume could 
cause cracks and gaps which would alter the effective thermal con¬ 
ductivity of the particle. Shrinkage and mechanical stress parame¬ 
ters are poorly known and more work is needed before they can be 
quantitatively incorporated into torrefaction particle simulations. 

Fig. 2b and d show the intraparticle temperature profiles at two 
radial positions and three times (t = 10, 20, 40 min) for the cases of 
final reactor temperatures of 235 and 289 °C, respectively. 
Intraparticle temperature gradients exists because of the internal 
heat transfer limitations which are expected to be significant as 
the Biot number (Bi = hL/k) is of order 1. For the Tf mat = 235 °C case, 
a positive temperature gradient across the particle develops for 
20 min and reaches a maximum during the drying stage at which 
point the surface of the particle is 50 degrees higher than the cen¬ 
ter. The experimental data shows that the temperature gradient 
remains positive throughout the entire hour. In the model, the cen¬ 
terline temperature very slightly exceeds the surface by approxi¬ 
mately 4 °C after 40 min. Despite this apparent qualitative 
difference, the difference between the model prediction and exper¬ 
imental centerline temperature is less than 8 °C which is well 
within the experimental error (±15 °C). 

For the T final =289°C case, a positive temperature gradient 
across the particle develops for 23 min where the surface of the 
particle is up to 85 degrees higher than the centerline. This positive 
gradient vanishes by 32.5 min and eventually becomes signifi¬ 
cantly negative. In fact, the centerline becomes hotter than the sur¬ 
face by up to 26 °C agreeing with centerline-surface temperature 
overshoot of up to 20 °C in the experiment. The cause for this phe¬ 
nomenon is the heat released from the exothermic torrefaction 
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Final reactor temperature (°C) 



Final reactor temperature (°C) 


Fig. 3. Effect of final reactor temperature on centerline temperature overshoot (a) 
and final (t = 60 min) mass/energy yields (b) for a 14 mm radius beech cylinder 
initially at 20.5 °C exposed to 10°C/min reactor heating rate, and a 27 °C initial 
reactor temperature. Experimental data in (a) from Van der Stelt [8], 


reactions. According to Fourier’s law, the heat generated at the 
inner layers can only conduct towards the outer layers of the par¬ 
ticle when a negative temperature gradient exists (dT/dr < 0). It is 
only when a sufficiently negative temperature gradient has devel¬ 
oped that a maximum in the centerline temperature is reached 
wherein the heat transfer by conduction and convection balances 
the heat released by reactions. During this period of thermal over¬ 
shoot, the heat transfer by convection of hot volatiles outwards 
through the particle is small - but significant - and contributes 
to heating the outer layers. 

4.1.3. Effect of temperature 

Van der Stelt repeated the experiments with the 28 mm diam¬ 
eter beech particles at final reactor temperatures ranging between 
200 and 280°C. In Fig. 3a, model predictions for the centerline tem¬ 
perature overshoot are compared with the experimental results. 
The overshoot is defined by the largest positive temperature differ¬ 
ence between the centerline and the reactor. At 200 °C the rate of 
torrefaction reactions is negligible producing few volatiles and 
therefore a very small (<1 °C) thermal overshoot. The magnitude 
of the overshoot rises rapidly with increasing temperatures and 
is certainly a non-linear effect. This is because the heat release 
and temperature are closely coupled: higher temperatures result 
in faster heat release increasing the temperature further. The over¬ 
shoot stabilizes itself when a sufficiently negative temperature 
gradient exists in the particle, thereby enabling conduction of the 
heat outward. 

Good quantitative agreement on the magnitude of the thermal 
overshoot as a function of temperature is reached between the 
model and experiment. Using a simplified particle model, Van 


der Stelt concluded that the temperature overshoot was propor¬ 
tional to the average volumetric heat release rate. Based on the 
thermal overshoot observed in the cylindrical particles, he 
estimated the cumulative heat release to be less than 250 kj/kg. 
The present modeling results corroborate this conclusion; the pre¬ 
dicted cumulative heat release is less than 180 kj/kg. As originally 
concluded by Turner et al. [16] the particle temperature profile 
during this overshoot period is extremely sensitive to the assumed 
thermochemistry model. In this case, the heat of reaction was not 
adjusted to fit the data. Instead its value depends solely on the 
composition of volatile and solid products which were fitted 
separately in [57] and whose enthalpies were estimated in [58], 
Fig. 3b shows the model predictions for the final (t = 60min) 
mass/energy yield for single 28 mm diameter beech wood particles 
as a function of reactor temperatures between 200 and 280 °C. For 
a given residence time and particle size, mass and energy yields 
decrease with temperature. Higher temperatures result in more 
rapid mass loss. The highly oxygenated volatiles released during 
torrefaction contain less energy than the original solid product 
[58], As a result the energy yield of the remaining solid product 
is always higher than the mass yield. This results in a product with 
greater energy density per unit mass. Van der Stelt did not measure 
the energy yield of the particles after the experiment so a compar¬ 
ison on this quantity is not available. Having thoroughly described 
and explained the phenomena associated with coupled thermal- 
kinetic effects at variable temperatures for relatively large cylindri¬ 
cal particles (2.8 cm diameter), the next section examines and val¬ 
idates the model predictions for various particle sizes. 


4.2. Comparison with the experiments of Basu and coworkers 


4.2.1. Description of experiments 

Basu and coworkers [22] performed torrefaction experiments 
using large cylindrical poplar (Populus deltoides) and particles with 
diameters between 0.476-2.54 cm and lengths of 0.8-6.5 cm and 
initial moisture contents of 6.39% w.b. The particles were sus¬ 
pended by a thermocouple from an electronic balance which mon¬ 
itored the weight of the sample. The thermocouple continuously 
recorded the temperature at the center of the particle (r = 0 cm). 
The convective reactor was a 42 mm diameter Quartz Wool Matrix 
reactor which was externally electrically heated. Nitrogen was 
used as an inert gas flowing at a rate of 1.2 L/min. In each case 
the higher heating value of the torrefied particle after the experi¬ 
ment was measured using a bomb calorimeter in order to evaluate 
the energy yield: 




m SJ HHV ss 

m s ,oHHVs, 0 ’ 


(19) 


where HHV S f is the final higher heating value of the torrefied parti¬ 
cle. The mass and energy yields were only reported at the end of the 
experiment at 60 min and the effect residence time was not studied. 


4.2.2. Effect of particle size 

In Fig. 4a, the model predictions for the centerline-reactor tem¬ 
perature overshoot as a function of particle size are compared with 
the experimental results from Basu et al. [22] and Dhungana [64], 
In all cases, the reactor temperature was held constant at 250 °C 
while the initial particle temperature was approximately 25 °C. 
The model predicts a 12.75 mm radius particle would experience 
an overshoot of 11.1 degrees which lies almost exactly between 
the two experimental measurements (7 and 18 °C). Both data sets 
(Basu et al. and Dhungana) refer to the same set of experiments 
done by the same group but mentioned in different publications. 
It is presumed that this discrepancy is because of the uncertainty 
on the actual reactor and particle centerline temperatures. 
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♦ Dhungana, 2011 
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Fig. 4. Effect of the particle size on the centerline-reactor temperature overshoot (a) particle size versus heat-up time (b) mass/energy yields 
Experimental data from Basu et al. [22], reactor temperature = 250 °C. 


particle size (c). 


Nevertheless, the trend is clearly visible in both model and 
experimental results: the magnitude of the centerline-reactor ther¬ 
mal overshoot increases quadratically with particle size. The theo¬ 
retical explanation for this trend was first analyzed by Van der Stelt 
[8], where a simplified quasi-steady state analysis of the energy 
conservation Eq. (12) showed that the maximum temperature dif¬ 
ference between the centerline and the surface of the particle 
should be proportional to the square of the particle radius (R in 
m 2 ), the volumetric heat release rate (-Q., in W/m 3 ) and inversely 
proportional to the effective conductivity (k in W m 1 K ') that is, 
AT C _ S oc -Q r R 2 /k. 

Fig. 4b shows the time for the centerline of the particle to heat 
up from its initial temperature to 200 °C. This heatup time is 
defined this way because torrefaction reactions do not initiate until 
200 °C. The model predicts the heatup time within 45 s of the 
experimental values. For the largest diameter particle, the heatup 
period takes more than a quarter (16.5 min) of the total residence 
time of the particle (60 min). 

The length of the heatup period is governed by two factors. The 
first is the rate of heat transfer to the particle which can be quan¬ 
tified by the characteristic time (in seconds) associated with solid 
particle heating. These scale linearly, pc p L/h, or quadratically, pc p . 
L 2 lk with size, depending on whether the heat transfer is limited 
by external or internal heat transfer rates, respectively. In this case, 
the Bi ranges from 0.5 to 5 for the smallest particle to largest par¬ 
ticle so that neither external nor heat transfer is rate limiting and 


the heatup time should be expected to scale with particle size with 
an order between 1 and 2. The second controlling factor during this 
period is the moisture content of the particle. It has been previ¬ 
ously been observed that the length of time for drying varies line¬ 
arly with the initial moisture content of moist wood cylinders 
subjected to convective heating [34], This occurs for several rea¬ 
sons: (1) liquid water has a more than double the specific heat 
capacity of biomass (4800 and -2300 J/kg K, respectively), which 
is not offset by its contribution to increasing the effective thermal 
conductivity [34], This lengthens the aforementioned characteris¬ 
tic times scales for heating. (2) Drying slows the particle heating 
rate because it is a strongly endothermic process and acts as a ther¬ 
mal sink. At certain locations in the particle, the thermal load for 
evaporation balances the rate of heat conduction inwards, and 
the temperature does not continue to rise until all remaining mois¬ 
ture has been evaporated. (3) The rapid outward convection of 
vapor from inner layers contributes to cooling the outer layers. 
Thus the ability of the model to predict heatup time well suggests 
that these key processes are described satisfactorily. 

Fig. 4c shows the model predictions for the mass/energy yields 
for the poplar particles of various sizes at 250 °C. Experimental 
results from Basu et al. demonstrate that the mass and energy 
yields decrease slightly with increasing particle size though the 
significance of this trend is not clear [65], Likewise the model pre¬ 
dicts no significant variation in the mass and energy yield with 
particle size. The energy yield is predicted to within 3% of the 
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Time (minutes) 


Fig. 5. Centerline temperature versus time (minutes) left axis for various particle sizes (radius = 2.38, 5, 7.5, 12.75 mm). Mass yield 
temperature = 250, 280 °C for (a) and (b) respectively. 


(right axis), reactor 


experimental energy yield for all but the smallest particle 
(r = 2.38 mm). The discrepancy could be attributed to the fact that 
the mass loss kinetics and thermochemical submodel are based on 
data from willow while in the experiment poplar was used. Unfor¬ 
tunately, the thermochemical data for poplar is lacking so the 
parameters for willow (a similarly deciduous tree) were used 
instead. 

The trends shown by the model and experiment cannot be 
explained by mere intuition. It would be expected that for a given 
residence time, larger particles would display less conversion 
(higher mass yield) than smaller particles. A larger particle 
spends more time heating up therefore experiences a shorter time 
at the reaction temperature, and should exhibit lower conversion 
(higher mass yield). The actual trends can only be explained by 
examining the predicted conversion profile of the different sized 
particles. In Fig. 5a and b the centerline temperature profile and 
mass yield are shown for reactor temperatures of 250 and 
280 °C, respectively. Initially, during the heat up period, the larger 
particles do in fact take longer to heat up and dry and therefore 
display lower centerline temperatures, and higher mass yields. 
The centerline of the 12.75 mm radius particle takes 30 min to 
reach the 250 °C reaction temperature while the centerline of 
the 2.38 mm radius particle takes only 6 min. During this time, 
the mass yield shows a strong dependence on particle size. For 
example at 20 min, the mass yield for the largest (12.75 mm 
radius) and smallest (2.38 mm radius) particles are 0.905 and 
0.985, respectively. However, once the centerline of the particle 
significantly exceeds 200 °C (the temperature at which torrefac- 


tion reactions initiate), the effect of the exothermic reactions 
becomes the controlling factor. Since the centerline temperature 
overshoot increases with particle size, larger particles exhibit 
higher temperatures and accelerated reaction rates. In the case 
of 250 °C, the reactions are accelerated enough so that the differ¬ 
ence in mass yield between the large and small particles 
decreases in time. By t = 60 min, the difference in conversion 
between particles of different sizes is very small. Once the ther¬ 
mal overshoot vanishes, the process no longer depends on the 
particle size (a relatively uniform temperature exists across the 
particle), and is controlled purely by kinetics. 

In the case of 280 °C (Fig. 5b), the exothermic temperature over¬ 
shoot effects are even larger and although the reactions are accel¬ 
erated strongly by this effect, the large particles still exhibit higher 
mass yield than the small particles throughout the process. At the 
end, t = 60min, the largest particle (12.75 mm radius) exhibits a 
mass yield nearly the same (0.62) as the smallest size particle 
(2.38 mm radius) which shows a mass yield of 0.615. In all cases, 
the behavior of the temperature and conversion profiles can be 
classified into three different periods: In the first period, the parti¬ 
cle heat-up is important and therefore the profiles are governed by 
the moisture content and heat transfer characteristics of the parti¬ 
cle and reactor system. In the second period, the exothermic over¬ 
shoot effects - which are sensitive to temperature, kinetics, and 
particle size - become the controlling factor. In the third period, 
the particle is at uniform temperature and the conversion profile 
is governed by the mass loss kinetics. The results highlight the 
importance of coupled, transient models which can account for 
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Fig. 6. (a) Heatup time versus wet basis moisture content and (b) dry mass yield versus time for two different particle sizes (R = 2.38,12.75 mm), and three different initial 
moisture contents (5,15, 30% w.b.). 


all three controlling phenomena (heat transfer, thermochemistry, 
and kinetics). 


4.3. Effect of moisture content 

In order to understand the effects of moisture content in more 
detail, a parametric modeling study was completed. The simulation 
conditions were chosen as representative of those utilized by Basu 
and co-workers (i.e. initially ambient temperature cylindrical pop¬ 
lar particles dropped into an isothermal reactor at 280 °C). Two 
variables were compared parametrically: the particle radius was 
varied between 2.38 and 12.75 (mm) and the initial moisture con¬ 
tent from 5 to 30 (% w.b.). In Fig. 6a, the heatup time (the time for 
centerline to reach 200 °C) is plotted versus the full range of vari¬ 
ables simulated). A wide range in the heatup times from as low as 
1.2 min for the smallest and driest particle (2.38 mm, 5% w.b) to 
35.4 min for the largest and wettest (12.75 mm, 30% w.b) particle 
is predicted by the model. The heating time varies linearly with 
moisture content and quadratically with particle size which agrees 
with results from the experimental drying results Di Blasi et al. 
[34], Therefore, the larger the particle the more sensitive the length 
of the drying time is to the initial moisture content. For example, a 
10% increase in moisture content lengthens the heatup period for a 
12.75 mm radius particle by 8 min while a 5 mm radius particle 
only takes 2.2 additional minutes for the same increase in moisture 
content. This sensitivity arises because the drying rate is largely 
governed by the temperature of the particle which - as discussed 
previously - depends on the particle size. Overall, the results sug¬ 
gests that when torrefying large (r > 10 mm) moist particles, atten¬ 
tion must be paid to the variations in moisture content between 
particles which will cause significant variations in conversion even 
for particles of the exact same size. 


In order to demonstrate these sensitivities more clearly, the dry 
mass yield profiles for two selected particle sizes of 2.38 mm and 
12.75 mm and moisture contents of 5, 15, 30% (w.b.) are shown 
in Fig. 6b. For all particles, there is an initial period where the 
dry mass yield is unchanged (i.e. torrefaction reactions have yet 
to initiate). This delay in the initiation of conversion is attributable 
to the difference in time needed to heat and dry the particles. For 
example, the large (12.75 mm) and wet (30%) particle takes more 
than 30 min before significant mass loss occurs. 

The conversion profile of the 2.38 mm particle is largely insen¬ 
sitive to the increased levels of moisture content. For this particle 
size, throughout the entire 60 min residence time, at no point does 
the difference in conversion level (quantified by the mass yield) 
between the dry (5% w.b.) and the moist (30% w.b.) vary by more 
than 7%. However, the conversion profiles of the large particle 
(12.75 mm) are extremely sensitive to increased moisture content. 
For example, at 35 min the dry (5% w.b.) particle is almost com¬ 
pletely torrefied (the mass yield is <70%) whereas the moist (30% 
w.b.) particle has yet to undergo significant torrefaction. This fur¬ 
ther illustrates that grinding of particles prior to torrefaction has 
a significant advantage: the process becomes much less sensitive 
to the initial moisture content. Although the temperature profiles 
are not shown, the magnitude of the exothermal overshoot is unaf¬ 
fected by moisture content. This is because the exothermal over¬ 
shoot does not initiate until the particle is completely dry. 

4.4. Implications for reactor design and process control 

Several important and novel conclusions regarding the effects 
of particle size on torrefaction have been elucidated by the present 
modeling results. Initially, the larger particles take longer than 
smaller particles to heat up to the reactor temperature. The 
disparity in time for heat up between small and large particles is 
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exacerbated by increasing moisture content. During this period, 
conversion decreases strongly with increasing particle size and 
increasing moisture content. This suggests that it may be advanta¬ 
geous to have a separate grinding and/or drying step prior to the 
torrefaction process which would help reduce the sensitivity of 
the process to particle size. During the second phase, the subse¬ 
quent effect of temperature overshoot caused by exothermic reac¬ 
tions accelerates the conversion of larger particles partially 
cancelling the effect of the delay caused by their slow heat up. 
Although the time at which the exothermal overshoot occurs is 
delayed by increasing moisture content, the magnitude of the over¬ 
shoot is unaffected by initial moisture content. As a result, the 
overshoot can be reduced to a minimum only by using a low torre¬ 
faction temperature 200-240 °C or small particle size. During the 
last thermal equilibrium stage of torrefaction, further conversion 
is governed solely by mass loss kinetics. 

The results demonstrate that the torrefaction of small and dry 
particles will be preferable to large and moist particles because 
of the formers’ more rapid and uniform conversion properties. 
However several economic and energetic trade-offs exists between 
performing drying and/or grinding prior to the torrefaction pro¬ 
cess. For example, one could propose a system configuration in 
which a grinding step to small ~2 mm particles occurs before a 
combined drying and torrefaction process. Although it shifts the 
thermal drying load to the torrefaction reactor, it negates the need 
for a separate drying step while also reducing the residence time, 
size, and cost of the reactors involved. The feasibility of this config¬ 
uration will ultimately depend on grinding equipment limitations 
as well as the initial moisture content of the raw biomass: mois¬ 
ture has been shown to increase the energy requirement and 
reduce throughput during grinding [5], For raw biomass types that 
can be delivered relatively dry and are easily comminuted - such 
as some herbaceous biomasses - this may be a suitable 
configuration. 

For feedstocks where this is not the case - such as woody bio¬ 
mass - the solution may be to design the torrefaction process so 
that the effects of particle size and moisture content on the overall 
conversion are minimized. This is done through two ways. One 
route is to meticulously characterize the feedstock to achieve 
incoming streams of uniform size and moisture content. Once 
known, the larger particles with higher moisture content could 
be treated separately from those with smaller sizes and lower 
moisture content. The other route is to modify the reactor heating 
conditions either through low temperatures and long residence 
times or through non-linear reactor heating profiles, reducing the 
sensitivity of the process to particle size. 

Overall when identifying the optimal configurations and resi¬ 
dence times, three competing factors must be quantified during 
torrefaction: the targets include (i) minimization of feedstock pre¬ 
treatment (i.e. grinding, sieving, drying, characterization) costs, (ii) 
minimization of reactor complexity, cost and residence time, and 
(iii) maximization of the energy yield of the solid products and 
overall efficiency. Using the coupled model described here to iden¬ 
tify optimal conditions which account for these three competing 
considerations will be the focus of future work. 


5. Conclusion 

A one-dimensional, coupled model for the torrefaction of single 
woody biomass particles has been presented and validated against 
published experimental data. The model is able to predict the ther¬ 
mal behavior as well as the consequent mass and energy yield of 
initially moist single particles undergoing torrefaction. The results 
show that the moisture content, particle size, residence time, and 
reactor temperature have strong and -in some cases nonlinear- 


effects on the mass and energy yield of the particle. The conversion 
profiles are characterized by three consecutive stages including (1) 
particle heat-up and drying, (2) exothermal overshoot, and (3) 
thermal equilibrium. The length of time required for particle 
heat-up and drying is determined by the external/internal heat 
transfer limitations and initial moisture content. As a result, con¬ 
version decreases with increasing particle size and moisture con¬ 
tent. During the second exothermal overshoot stage, heat 
released from chemical reactions causes significant intraparticle 
gradients and also accelerates the conversion of large particles. 
While the assumed kinetic and thermochemical parameters are 
clearly crucial to determining the magnitude of this exothermal 
overshoot stage, these were not adjusted or tuned during the val¬ 
idation of this particle model. Lastly, when thermal equilibrium 
is reached, the particle conversion is governed by solely by mass 
loss kinetics. The validated model can serve as a useful tool in 
future work to assess the reactor conditions (temperature, resi¬ 
dence time) and feed particle sizes necessary to achieve a homoge¬ 
neous and sufficient degree of conversion as well as for quantifying 
the trade-offs between various aspects of the torrefaction process. 


Acknowledgements 

The authors gratefully acknowledge BP for funding this 
research. The authors would also like to thank Dr. Michiel van 
der Stelt for sending a copy of his dissertation. 


References 

[1] Bergman PCA, Boersma AR, Kiel JHA, Zwart RWR. Torrefaction for biomass co¬ 
firing in existing coal-fired power stations. Energy Res Center Neth 2005. 

[2] Arias B, Pevida C, Fermoso J, Plaza MG, Rubiera F, Pis JJ. Influence of torrefaction 
on the grindability and reactivity of woody biomass. Fuel Process Technol 
2008;89:169-75. 

[3] Prins MJ. Thermodynamic analysis of biomass gasification and torrefaction. 
PhD thesis, Eindhoven University of Technology; 2005. 

[4] Hakkou M, Petrissans M, Gerardin P, Zoulalian A. Investigations of the reasons 
for fungal durability of heat-treated beech wood. Polym Degrad Stab 
2006;91:393-7. 

[5] Miao Z, Grift TE, Hansen AC, Ting KC. Energy requirement for lignocellulosic 
feedstock densifications in relation to particle physical properties, preheating, 
and binding agents. Energy Fuels 2012;27:588-95. 

[6] Prins MJ, Ptasinski KJ, Janssen FJJG. Torrefaction of wood: Part 1. Weight loss 
kinetics. J Anal Appl Pyrol 2006;77:28-34. 

[7] PengJH, Bi HT, Sokhansanj S, Lim JC. A study of particle size effect on biomass 
torrefaction and densification. Energy Fuels 2012;26:3826-39. 

[8] Van der Stelt MJC. Chemistry and reaction kinetics of biowaste torrefaction. 
PhD thesis, Eindhoven University of Technology; 2011. 

[9] PengJH, Xiaotao TB, Lim J, Sokhansanj S. Development of torrefaction kinetics 
for British Columbia softwoods. Int J Chem React Eng 2012;10. 

[10] Chen W-H, Kuo P-C. Isothermal torrefaction kinetics of hemicellulose, 
cellulose, lignin and xylan using thermogravimetric analysis. Energy 2011; 
36:6451-60. 

[11] Perre P, Remond R, Turner I. A comprehensive dual-scale wood torrefaction 
model: application to the analysis of thermal run-away in industrial heat 
treatment processes. Int J Heat Mass Transfer 2013;64:838-49. 

[12] Chan W-CR, Kelbon M, Krieger BB. Modelling and experimental verification of 
physical and chemical processes during pyrolysis of a large biomass particle. 
Fuel 1985;64:1505-13. 

[13] Gronli MG, Melaaen MC. Mathematical model for wood pyrolysis: comparison 
of experimental measurements with model predictions. Energy Fuels 2000;14: 
791-800. 

[14] Lu H, Ip E, Scott J, Foster P, Vickers M, Baxter LL. Effects of particle shape and 
size on devolatilization of biomass particle. Fuel 2010;89:1156-68. 

[15] Haseli Y, van Oijen JA, de Goey LPH. A detailed one-dimensional model of 
combustion of a woody biomass particle. Bioresour Technol 2011;102: 
9772-82. 

[16] Turner 1, Rousset P, Remond R, Perre P. An experimental and theoretical 
investigation of the thermal treatment of wood ( Fagus sylvatica L.) in the range 
200-260 °C. Int J Heat Mass Transfer 2010;53:715-25. 

[17] Kadem S, Lachemet A, Younsi R, Kocaefe D. 3d-Transient modeling of heat and 
mass transfer during heat treatment of wood. Int Commun Heat Mass Transfer 
2011;38:717-22. 

[18] Patuzzi F. Gasparella A, Baratieri M. Thermochemical and fluid dynamic model 
of a bench-scale torrefaction reactor. Waste Biomass Valorization 2013:1-9. 



















